Atomic Carbon and CO Isotope Emission in the Vicinity of DR15 
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ABSTRACT 

We present observations of the 3 P\- 3 Pq fine structure transition of atomic carbon [CI], the 
J=3-2 transition of CO, as well as of the J=l— transitions of 13 CO and C 18 toward DR15, an 
HII region associated with two mid- infrared dark clouds (IRDCs). The 13 CO and C 18 J=l-0 
emissions closely follow the dark patches seen in optical wavelength, showing two self-gravitating 
molecular cores with masses of 2000 Mq and 900 Mq, respectively, at the positions of the 
catalogued IRDCs. 

Our data show a rough spatial correlation between [CI] and 13 CO J=l— 0. Bright [CI] emission 
occurs in relatively cold gas behind the molecular cores, neither in highly excited gas traced by 
CO J=3-2 emission nor in HII region/molecular cloud interface. These results are inconsistent 
with those predicted by standard photodissociation region (PDR) models, suggesting an origin 
for interstellar atomic carbon unrelated to photodissociation processes. 

Subject headings: HII regions — ISM: individual (DR 15) — ISM: atoms — ISM: molecules 
1. INTRODUCTION 
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The amount and distribution of atomic carbon 
is important for understanding of the physical and 
chemical structure, formation precesses, and ther- 
mal balance of interstellar molecular clouds. Since 
the ionization potential of neutral atomic carbon 
(11.266 eV) is close to the dissociation energy of 
CO (11.09 eV), it (C°) is thought to be confined in 
a thin layer near the surfaces of molecular clouds. 
Chemical models of molecular clouds exposed to 
ultraviolet (UV) radiation demonstrate the con- 
centration of C° in the range Ay =3-6 mag (Tie- 
lens & Hollenbach 1985; van Dishoeck & Black 
1988; Hollenbach, Takahashi, & Tielens 1991). 

The 3 Pi- 3 Po fine-structure line of neutral 
atomic carbon [CI] (492.160651 GHz, Yamamoto 
& Saito 1991) has been known to be a good tracer 
of interstellar atomic carbon. After the first de- 
tection (Phillips et al. 1980), the rapid advance 
in receiver sensitivity in the submillimeter wave- 
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length has resulted in the 3 Pi- 3 Pq line now being 
well studied in many Galactic sources. Large-scale 
observations have revealed that [CI] emission is 
widespread throughout molecular clouds (Plume, 
Jaffc, & Keene 1994; Plume et al. 1999; Maezawa 
ct al. 1999; Ikeda et al. 1999), being not confined 
in vicinities of ultraviolet sources. Observations on 
scales of tens of degrees has shown that [CI] emis- 
sion is widespread also on galactic scales (Wright 
et al. 1991; Bennett et al. 1994). 

There have been some arguments (Keene et al. 
1985; Schilke et al. 1993; Schilke et al. 1995; 
Maezawa et al. 1999; Ikeda et al. 1999) that 
C° exists also within deep inside molecular clouds 
where UV radiation is attenuated. Models of time- 
dependent chemistry (Suzuki 1979; Leung, Herbst, 
& Huebner 1984; Suzuki et al. 1992; Lee et al. 
1996) predict that C° can be abundant in the in- 
terior of a newly formed dense core, since the con- 
version timescale of C° to CO is comparable to the 
dynamical timescale of the core. Even in steady 
states, some models predict that, under certain 
circumstances, atomic carbon can be abundant in 
the absence of UV radiation (Pineau des Forets, 
Roueff, & Flower 1992; Le Bourlot, Pineau des 
Forets, & Roueff 1993; Le Bourlot, Pineau des 
Forets, Roueff, & Flower 1995; Le Bourlot, Pineau 
des Forets, & Roueff 1995; Lee, Roueff, Pineau des 
Forets, Shalabiea, Terzieva, & Herbst 1998). Such 
models give two stable solutions of chemical equi- 
librium, the low ionization phase (LIP) and the 
high ionization phase (HIP), at low densities. 

[CI] line observations from objects with high 
visual extinction (Ay) may be crucial to answer 
the question whether C° really exists deep inside 
molecular clouds. In this study, we chose DR15, 
an HII region in the Cygnus-X region, as a tar- 
get. DR 15 is associated with two mid-infrared 
dark clouds (IRDCs; Egan et al. 1998), which be- 
long to a population of compact objects with sig- 
nificant mid-infrared extinctions (A^^m > 2 which 
corresponds to Ay ;> 200). H2CO observations 
(Carey et al. 1998) have shown that IRDCs are 
cold (T k < 20 K) and dense [n(H 2 ) > 10 5 cm" 3 ]. 
Lack of mid- to far-infrared emission associated 
with these clouds suggests that they do not con- 
tain currently forming high-mass stars. DR15 it- 
self is formed by one or two BO ZAMS star(s), 
which is observed as the far-infrared source FIR 1 
(Odcnwald et al. 1990), having a total luminosity 



of ~ 3 x 10 6 Lq . Assuming that all the far- infrared 
flux of FIR 1 was initially emitted in the ultravio- 
let (6-13.6 eV), an UV enhancement factor Go is 
estimated to be ~ 550 at 1 pc from the ionizing 
star(s). 

This paper presents the results of CI 3 Pi- 3 Pq, 
CO J=3-2, 13 CO J=l-0 and C 18 J=l-0 obser- 
vations of the DR 15 region. These data are used 
to study the C° column density and C°/CO abun- 
dance ratio, and their spatial variations, especially 
in relation to the UV source and the IRDCs. 

2. OBSERVATIONS 

2.1. Mt. Fuji Observations 

The CI 3 Pi- 3 P (492.160651 GHz) and CO 
J=3-2 (345.795991 GHz) observations were per- 
formed with the Mount Fuji submillimeter-wave 
telescope (Sckimoto et al. 2000) in 1999 March 
and December, and 2000 January. The diameter 
of the telescope is 1.2 m, which gives half-power 
beamwidths of 2'. 2 at 492 GHz and 3'.1 at 346 
GHz. The pointing of the telescope was checked 
and corrected by observing the sun and the full 
moon every month, and its accuracy was main- 
tained within 20"(rms) during the observing sea- 
son. The beam efficiencies including radome loss 
(ilmoon) were measured by observing the full moon 
to be 72 % (492 GHz) and 75 % (346 GHz), and 
the main-beam efficiencies are estimated to be 45 
% (492 GHz) and 67 % (346 GHz). The telescope 
efficiencies and beamwidths at the observed fre- 
quencies are summarized in Table 1. 

The telescope is equipped with three SIS re- 
ceivers (346 GHz/ 492 GHz/ 809 GHz) operated 
in the double-sideband (DSB) mode. The typical 
DSB system noise temperatures during the obser- 
vations were 700 K for 492 GHz and 500 K for 
346 GHz. Calibration of the antenna temperature 
was accomplished by chopping between an ambi- 
ent temperature load and the sky. Reproducibility 
of intensities were checked by monitoring the cal- 
ibration sources, Orion KL and DR21, and were 
found to be reliable within 9 % (la) at 492 GHz 
and 11 % (la) at 346 GHz 1 during the DR15 
observations. All spectra were obtained with a 
1024 channel acousto-optical spectrometer (AOS) 
which covers an instantaneous bandwidth of 700 
MHz with a spectral resolution of 1.6 MHz. These 
correspond to a 430 km s _1 velocity coverage and 
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an 1 km s 1 velocity resolution at the [CI] line 
frequency. 

The [CI] observations were performed in a 
frequency-switching mode. Switching frequency 
was set to 80 MHz and intervals to 10 seconds. 
The on-source integration times was 100 s for 
each position and rms noise was 0.1 K in T\. 
The CO J=3-2 observations were performed by 
position switching to a clean reference position, 
(01950A950) = (20 h 16 m 00 s ,+43°00'00"). Typ- 
ical on-source integration time was 350 s which 
yields rms noise of 0.2 K in T'\. We mapped a 
27'x27' area in the [CI] line, and a 42'x30' area in 
the CO J=3-2 line, on a 1'.5 grid. We subtracted 
baselines of the spectra by fitting fourth-order 
polynomial lines. 

2.2. NRO 45 m Observations 

The 13 CO and C 18 J=l-0 observations have 
been made with the NRO 45 m telescope in Jan- 
uary 2000. The NRO 45 m telescope has a 17"±1" 
full width at half maximum beam at 110 GHz. 
Pointing errors were corrected every two hours by 
observing the SiO maser source UX-Cyg with the 
HEMT amplifier receiver at 43 GHz. The point- 
ing accuracy of the telescope was good to < 3" in 
both azimuth and elevation. 

We used the 2x2 focal-plane array SIS receiver 
S115Q. Typical system noise temperature was 400 
K during the observations. Calibration of the 
antenna temperature was accomplished by a stan- 
dard chopper-wheel method. Since the mixers in 
S115Q work in a quasi-single sideband mode, we 
scaled the antenna temperature for each channel 
referring to that taken with the single-beam SIS 
receiver S100 equipped with a quasi-optical im- 
age rejection filter. Every observing day, we took 
spectra of DR21(CO) (Richardson et al. 1986), 
at (ai 9 5o,<5i95o) = (20 h 37 m 09.0 s , +42°09'00") 
near its transit with each channel of S115Q, 
and calculated scaling factors by comparing with 
TJ( 13 CO) = 12.8 K or TJ(C ls O) = 2.12 K . We 
used the high-resolution acousto-optical spectrom- 
eters (AOS-Hs) each of which covers an instanta- 
neous bandwidth of 40 MHz with a spectral reso- 



Intensity reproducibility of Mt. Fuji submillimcter-wave 
telescope is limited by the effect of optical depth and vari- 
ation in the temperature difference between the chopper 
inside the radome and the outside atmosphere. 



lution of 37 kHz. At the 13 CO J=l-0 frequency, 
these correspond to a 110 km s^ 1 velocity coverage 
and a 0.1 km s _1 velocity resolution, respectively. 

All data were obtained by position switch- 
ing to a reference position, (0:1950, £1950) = 
(20 h 45 m 13 s ,+42°38'52"). To conduct the sur- 
vey effectively, three on-source positions were ob- 
served for one reference position observation. We 
mapped a 20' x 20' area in the 13 CO line, and more 
restricted area in the C 18 line, on a 34" grid. We 
also made a 17" grid 13 CO mapping in 2'. 5 x 2'. 5 
areas around the two mid-IR dark clouds and the 
CO J=3-2 emission peak. Typically 20 seconds 
on source integration gave spectra with rms noise 
of 0.5 K (TjJ) at a velocity resolution of 0.1 km 
s _1 . We subtracted baselines of the spectra by 
fitting first or third-order polynomial lines. 

3. RESULTS 

The integrated intensity map in the [CI] and 
CO J=3-2 line are shown in Fig.l with the optical 
image. Figure 3 shows 13 CO J=l-0 and C 18 
J=l-0 integrated intensity maps. All the [CI] and 
CO maps are presented in units of T\. In this 
paper, we adopt 1.0 kpc as the distance to DR15 
(Carey et al. 1998). At this distance, the beam 
sizes of the Mt. Fuji submillimeter-wave telescope 
correspond to 0.64 pc (492 GHz) and 0.90 pc (346 
GHz), respectively. 

3.1. CI Distribution 

The [CI] emission is extended over the map with 
patchy morphology being not restricted to a region 
adjacent to the HII region, DR15. It roughly fol- 
lows dark patches in the optical image. An intense 
[CI] lump with a size of 8'x6' (2.3x1.7 pc 2 ) lies far- 
ther from the ionization front than the two IRDCs 
do. The [CI] lump has a 'hump' in its southern 
portion, near the position of G79.34+0.33, while 
G79. 27+0.38 has no [CI] counterpart. The [CI] 
lump has an weak extension to the south which 
follows the brightest part of the optical nebula. 

3.2. CO J=3 2 Distribution 

The CO J=3-2 emission is also extended over 
the map with more smooth distribution than [CI] . 
It strongly peaks at the position of DR15 FIR 1 
(Odcnwald et al. 1990), which may coincides with 
the ionizing star(s), and exhibits a poor spatial 
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correlation with the [CI] distribution. The dis- 
placement of the CO J=3-2 peak to that of [CI] is 
about 7' (2 pc). An western, weak intensity emis- 
sion hump may correspond to G79.27+0.38, while 
G79.34+0.33 is included in the bulk emission as- 
sociated with DR15. 

3.3. 13 CO and C ls O Distributions 

The 13 CO J=l-0 and C 18 J=l-0 maps show 
two prominent cores with diameters of 2'-3' (0.6- 
0.9 pc), which spatially coincide with the IRDCs, 
G79.27+0.38 and G79.34+0.33. Both the 13 CO 
and C 18 emissions closely follow dark patches in 
the optical image. The most intense CO peak lies 
at the G79.34+0.33 core, being associated with 
a CI counterpart. A low intensity envelope of 
the 13 CO emission extends over the upper half of 
the map, showing a rough spatial correlation with 
the [CI] lump (Fig. 4). An weak 13 CO emission 
peak at (ai 950 , /3i 95 o) = (20 h 30 m 22 s , +40°13') cor- 
responds to the northwest corner of the [CI] lump. 
No CO counterpart of the [CI] peak was found. 
Two small clumps at (ai 950 , /3i 950 )^(20 h 30 m 40 s , 
+40°06') and another larger clump in the south 
correspond to the CO J=3-2 peak and its south- 
ern extension. Since the two small clumps coincide 
with a band of obscuration separating the optical 
emission nebula, they must be in front of the HII 
region. 

The two molecular cores which correspond to 
G79. 27+0.38 and G79.34+0.33 have LTE masses 
of 2.0 x 10 3 M and 9.1 x 10 2 M , respectively, 
if we assume T O x = 20 K (see §4.2). These masses 
are similar to their virial masses, 1.3 x 10 3 Mq 
and 9.2 x 10 2 Mq, suggesting that the cores are 
gravitationally bound. The mean densities derived 
from LTE masses are, n(H 2 ) ~ 1.2 x 10 5 cm~ 3 and 
~ 1.3 x 10 5 cm~ 3 , which are consistent with the 
densities determined from LVG analysis of H 2 CO 
lines (Carey et al. 1998). 

3.4. Integrated Line Intensities 

Figure 5 shows correlation plots between the 
line intensities integrated over the range Vlsr = 
—7.5 km s _1 to +7.5 km s _1 . All the data are 
smoothed with the FWHM=2' Gaussian weighting 
function, and resampled to the same regular grids 
with 1.5' spacings. There is a rough correlation 
between the [CI] intensity and the 13 CO J=l 



intensity (correlation coefficient, 7? = 0.62; Fig. 4a). 
A least-squares regression to a linear relation gave 
a following relation; 

J T MB (Cl)dV = (0.66 ± o.4 6 )+(0.25 ±0 .o2) J T MB ( 13 CO)dV. 

(1) 

No correlation between the [CI] intensity and 
the C 18 ,7=1-0 intensity has been found (7? = 
0.35; Fig. 46). The C 18 intensity is well corre- 
lated with the 13 CO intensity (7? = 0.85; Fig.4c). 
The CO J=3-2 intensity shows a poor correlation 
with the [CI] intensity (7? = 0.42; Fig.4d). This 
seems to be remarkable, because the both line 
emissions are generally thought to arise from UV- 
irradiated cloud surfaces (e.g., Warin, Benayoun, 
& Viala 1996). 

4. DISCUSSION 

4.1. Column Densities 

Since the critical densities of CI 3 Pi~ 3 Pq and 
CO J=l-0 lines are similar to each other (~ 10 3 
cm -3 ), we expect similar excitation temperatures 
(T ox ) for the both transitions if C° and CO reside 
in the same or adjacent layers. We have calcu- 
lated the C° and CO column densities by assum- 
ing homogeneous clouds in local thermodynamic 
equilibrium (LTE). 

The total beam-averaged atomic carbon column 
density is given by: 



N(C°) = 1.98 x 10 15 
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where Q(T CX ) is the ground state partition func- 
tion for the carbon atom 



Q(T ox ) = l + 3exp(- 
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E 1 is the energy of the J=l level (75i/fc = 23.6 K), 
and 7?2 is the energy of the 7=2 level (E\/k = 62.5 
K). The factor rci/(l — e~ TCI ) (> 1) is a correction 
factor for the velocity-averaged optical depth, and 
is unity at the optically thin limit (rci <1). We 
substitute the optical depth at the peak velocity 
for the velocity-averaged optical depth. This sub- 
stitution makes derived column density an upper 
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limit, while the optically thin assumption gives a 
lower limit. The peak CI optical depth is given 
by: 



tci = -In 



1 - 



T 



MB 



*?/(./„ (Tex) - J„(Tbb)) 



(4) 



where 77/ is the beam- filling factor, and J v is the 
radiation temperature 



UT) = 



hv/k 



exp(hu/kT) 



(5) 



We assumed r]f = l in the CI optical depth estima- 
tion. Arguments based on observations of the CI 
3 P 2 - 3 Pi/ 3 Pi- 3 P ratio suggest that the CI 3 J>i- 
3 Po opacity is usually < 1 (Zmuidzinas et al. 1988; 
Genzel et al. 1988) . Our observations also gave 
tci < 1, even for the [CI] intensity peak position 
(tqi = 0.6). The uncertainty in the estimated C° 
column density brought by the optical depth cor- 
rection may be less than 33 % . Although the 
sub-beam filling effect may raise the optical depth, 
it does not make a significant revision to the col- 
umn density, so far as the CI line is optically thin. 
The calibration error of submillimeter line inten- 
sity also brings an uncertainty in the column den- 
sity (9 %). 

The general equation for estimating the CO col- 
umn density is given by: 



N(CO) = 5 ; 5 il 1 i ° 1 2 7 / Tmb<1V T ex exp ( E 



(i/[GHz])" 



1 



J„(T] 



BBi 



J v (T c: 



kT cx 

TCO 



1 - e- r co 



where E u is the upper level energy, tco/(1 — 
e -Tco) i s a correction factor for the velocity- 
averaged optical depth of CO line (tco)- The CO 
optical depths are estimated by the 13 CO/C 18 
intensity ratios assuming the isotopic abundance 
ratio, [ 13 CO]/[C 18 0]= 7.4 (Langer & Penzias 
1990), and equal beam-filling factors. We used 
the C 18 integrated intensity to estimate the 
CO column density with the isotopic abundance 
[CO]/[C 18 0]=392 (Langer & Penzias 1990). The 
observed 13 CO/C 18 integrated intensity ratio 
(> 4) implies that the C 18 J=1~0 transition is 
usually optically thin (r<;0.2). Thus, uncertainty 
in the estimated CO column density brought by 
the C 18 optical depth correction is less than 12 



(6) 



% . The excitation temperature is estimated from 
the 13 CO J=l-0 peak antenna temperature using 
the relationship: 

t mb ( 13 co) = nj [J„(T ex ) - J V (T BB )} (1 - e - ri3 ) , 

(7) 

where 77/ is the beam-filling factor (assumed to be 
= 1), T13 is the 13 CO optical depth at the peak 
velocity. The sub-beam filling effect may not be 
not significant, especially for two molecular cores, 
since their mean LTE densities are consistent with 
the LVG densities (§3.3). 

Table 2 lists the C° and CO column densities 
toward four representative directions in DR15, the 
[CI] peak, the centers of IRDCs (G79. 37+0.33, 
G79. 27+0.38), and CO J=3-2 peak. The data are 
smoothed with the FWHM=2' Gaussian weighting 
function. Excitation temperatures are typically 20 
K except for the CO J=3-2 peak (T cx ~40 K). The 
beam- averaged C°/CO column density ratios to- 
ward IRDCs are ~ 0.02, which are much smaller 
than those observed in the cloud cores (0.16-0.32, 
Plume et al. 1999; Ikeda et al. 1999). The CO 
J=3-2 peak, the position coincides with DR15, 
also shows low value of N(C°)/N(CO) (~ 0.03). 
In the [CI] peak, we measure 7V(C°)/iV(CO) to be 
0.15, which is similar to the typical value for the 
bulk of giant molecular clouds and nearby dark 
clouds (0.10-0.21, Ikeda et al. 1999; 0.08-0.25, 
Maezawa et al. 1999), although the line-of-sight 
extinction is far larger than those of them. 

If we assume T cx = 18 K and optically thin, 
the total C° number in the 13 CO mapping area 
amounts to 1.8 xlO 55 atoms. This total C° num- 
ber varies by 25 % for the excitation temperature 
range from 15 K to 40 K. The optically thin as- 
sumption makes us underestimate the C° number 
by up to 33% . The total CO number in the same 
area estimated under the optically thin assump- 
tion is 2.8 x 10 56 molecules. This is an underesti- 
mation by up to 50 %, and varies by a factor of 2 
for 15 < T cx < 40 K. The C°/CO ratio for the bulk 
of the DR15 cloud is 0.064 ± 0.035. 

4.2. Origin of [CI] Emission 

4.2.1. Does UV from DR15 Dominate C° For- 
mation ? 

Low C°/CO ratios toward cloud cores are re- 
flected in the [CI] vs. 13 CO intensity correla- 
tion plot (Fig. 5 a). This trend has been found 
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in the data obtained by the previous observations 
(Tauber et al. 1995; Plume et al. 1999, Ikeda et 
al. 1999), and is qualitatively consistent with the 
steady-state PDR models, that C° resides only in 
thin surfaces of UV-irradiated molecular clouds. 
Poor correlation of [CI] intensity with C 18 J=l- 
intensity (Fig. 5b) also implies that [CI] emis- 
sion arises from cloud surfaces and C 18 emission 
prefers cloud cores. 

However, the [CI] intensity peaks farther into 
the molecular cloud from the ionization front be- 
hind the IRDCs (Fig. 16). The HII region / CO 
core / CI lump arrangement itself is inconsistent 
with the expectation from steady-state, plane- 
parallel PDR models. In M17 and S140, sim- 
ilar arrangements have been reported by Keene 
et al. (1985). We have found the C+'/CO/C 
arrangement also in the Orion KL region (Ikeda 
et al. 2001) and the p-Oph main cloud (L1688; 
Kamegai et al. 2001). If the [CI] peak is shad- 
owed by G79.37+0.33, the visual extinction be- 
tween the ionization front and the [CI] peak po- 
sition exceeds 80 mag. The [CI] peak direction 
itself has a line-of-sight extinction of 20-30 mag. 
On the other hand, models for plane-parallel pho- 
todissociation regions have shown that the layer 
with N(C°)/N(CO) >0.1 has depth smaller than 
7 mag. Clumpy structure may help UV pho- 
tons penetrate deep into molecular clouds. How- 
ever, if we assume n(H 2 ) = 10 4 ~ 5 cm~ 3 and an 
empirical relation N(H 2 ) = 1.6 x 10 21 A V cm~ 2 , 
N(C°)/N(CO) < 0.1 requires clumps with sizes 
larger than 0.4-4 pc, which is larger than the res- 
olution of our observations (NRO 45 m: 0.16 pc; 
Mt. Fuji: 0.64 pc). 

In addition, the integrated intensity of the CO 
J=3-2 line, which also arises from UV-irradiated 
cloud surfaces (e.g., Gierens, Stutzki, & Win- 
newisser 1992), shows poor spatial correlation 
with the [CI] intensity (Fig.l). The CO J=3-2 
intensity strongly peaks at DR15 FIR 1, having 
an extent of 5' x 10' (1.5 x 3 pc 2 ), while the [CI] 
intensity peaks farther into the molecular cloud 
from the ionization front than either of the J=l- 
line of the CO isotope does. The situation is 
well illustrated in Fig. 6 which shows a plot of in- 
tegrated intensities along a line passing through 
the [CI] peak and DR15A 2 . A hump at DR15 
in the [CI] intensity is the southern extension of 
the [CI] lump (§3.1), and it could correspond to 



the emission from thin photodissociated layer near 
surfaces of UV-irradiated clouds. 

The variation of the density could account for 
the observed arrangement. The lower density pro- 
duces a lower temperature photo-dissociated layer, 
which has low rates of OH production and, con- 
sequently, increases C° by shutting off the CO 
formation route through OH (e.g., Hollenbach, 
Takahashi, & Tielens 1991). Following this idea, 
the observed C° column densities require gas den- 
sity around DR15 at least an order of magnitude 
higher than the [CI] peak and the dense molecular 
core G79.37+0.33. G79. 27+0.38 has the C° col- 
umn a factor of 3 lower than that of G79.37+0.33, 
although the both have similar densities (§3.3) and 
UV environments. Apparently, G79.37+0.33 has 
a [CI] counterpart, while G79.27+0.38 does not. 
The mean H 2 density toward the CO 3-2 peak, 
which is estimated from the 13 CO column density, 
is not so high, ra(H 2 )<;2 x 10 4 cm~ 3 . These sug- 
gest that the density is unlikely a key parameter 
to determine the C° column density in the DR15 
region. 

The C° column density toward the CO 3-2 peak 
is estimated to be 9.3 xlO 16 cm~ 2 . According to 
Hollenbach, Takahashi, & Tielens (1991), a plane- 
parallel PDR with Go = 500 and n = 10 4 cm" 3 
has a C° column density of 2.5 x 10 17 cm~ 2 . This 
value differs from the observed value by more than 
a factor of 2. Indeed, Go may be higher than 500 at 
the position close to DR15A, the discrepancy be- 
tween theory and observation would be more sig- 
nificant. The hypothesis of unresolved, high den- 
sity clumps, does not change the situation, since 
it decreases the CI column density of each clumps 
while it increases the surfaces to be consistent with 
the observed 13 CO intensity. We therefore con- 
clude that the steady-state PDR models are not 
able to give reasonable interpretations to the ob- 
served CO 3-2(HII)/CO/CI arrangement and the 
deficit of C° toward DR15. 

4-2.2. Origin of Atomic Carbon Unrelated to 
Photodissociation 

The similarity in the distributions of the [CI] 
and 13 CO lines may be an indication that these 
lines arise from the same spatial regions. This 

2 A compact radio source located at the brightest part of FIR 
1 (Collcy 1980). 
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implies that C° basically coexist with molecular 
gas, except for the dense regions in molecular 
clouds, such as molecular cores G79.37+0.33 and 
G79.27+0.38. This idea is supported by the re- 
sults of our large-scale CI surveys of nearby molec- 
ular clouds (e.g., Ikeda et al. 1999). To account 
for the coexistence of C° and CO, some revisions 
to the steady-state PDR models or the application 
of models unrelated to photodissociation may be 
inevitable. It may be constructive to examine the 
other ideas for the origin of atomic carbon unre- 
lated to photodissociation processes. 

Timc-dcpendcnt chemistry models (e.g., Lee et 
al. 1996; Storzer, Stutzki, & Sternberg 1997; 
Suzuki et al. 1992) can explain the observed [CI] 
and CO distributions. The dynamical timescale of 
a dense core with n(H 2 ) = 10 5 cm -3 is td yn — 10 5 
years, while its chemical timescale (timescale for 
conversion from C° to CO) is t c h em <, 10 4 years. 
Thus, the carbon in a high-density core is incor- 
porated into CO completely before the core col- 
lapse. In the less dense region with n(H 2 ) = 10 4 
cm~ 3 , the dynamical timescale and the chemical 
timescale become comparable [(3^4)xl0 5 years], 
and hence C° can coexist with CO for a while even 
after UV radiation is shielded. The Cl-rich cloud 
found in Heilcs Cloud 2 (Maezawa et al. 1999) 
may be the typical of such a chemically young 
molecular cloud. It is likely that high-density 
cores, G79.37+0.33 and G79.27+0.38, are chem- 
ically evolved, while less dense envelope is chemi- 
cally young. 

Another idea is based on the chemical models 
with bistable solutions. These models predict that 
there exist two phases in the chemical equilibrium 
state at low densities (n.H<j5500 cm~ 3 ; Pineau des 
Forets, Roucff, & Flower 1992). The two stable 
solutions referred to as the high ionization phase 
(HIP) and the low ionization phase (LIP) lead 
into two different C°/CO ratios in clouds with 
identical density and temperature. This means 
that the gas-phase abundance of C° can be large 
(C°/CO=0.1-0.2) in HIP, even in the absence of 
a UV field. For the cloud cores, where the density 
exceeds the critical value for bistability, the sys- 
tem necessarily take a solution with a low C°/CO 
ratio (< 0.01). The chemical bistability can be 
applicable to the DR15 case, if we consider the 
[CI] peak as HIP and the cores of G79. 37+0.33 
and G79. 27+0.38 as LIP. The bistability model 



can also account for the ubiquitous value of the 
C°/CO ratio (0.1-0.2) over the bulk of molecular 
clouds, although its applicability is still controver- 
sial. 

5. CONCLUSIONS 

We have studied the distributions of the CI 
3 Pi- 3 P line, the CO J=3-2 line, as well as of 
the J=1~0 transition of 13 CO and C 18 toward 
DR15, an HII region associated with two IRDCs. 
Neutral atomic carbon is extended over the molec- 
ular clouds traced by the 13 CO J=l-0 line. The 
two dense molecular cores with masses of 2.0 x 10 3 
Mq and 9.1 x 10 2 Mq, which correspond to the 
IRDCs, have been found in the 13 CO and C 18 
data. The [CI] integrated intensity peaks far- 
ther into the molecular cloud from the ioniza- 
tion front behind the IRDCs. The column den- 
sity ratio N(C°)/N(CO) toward the [CI] peak is 
0.15, which is similar to those observed in giant 
molecular clouds and nearby dark clouds, while 
the dense cores (IRDCs) exhibit very low ratios 
(~0.02). The CO J=3-2 emission strongly peaks 
at the DR15 FIR 1, exhibiting a poor spatial cor- 
relation with the [CI] distribution. The CO J=3- 
2(HII)/CO/CI arrangement and the deficit of C° 
toward DR15 are inconsistent with the prediction 
of PDR models, and suggest an origin for inter- 
stellar atomic carbon unrelated to photodissoci- 
ation. Models with time-dependent chemistry or 
chemical bistability could account for the observed 
H+/CO/C arrangement and the C°/CO ratios as 
well. 

We are grateful to the NRO staff for excellent 
support in the 45 m observations. This study 
is supported by Grant-in-Aid from the Ministry 
of Education, Science, and Culture (07CE2002, 
11304010, and 12740119) 
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Fig. 1. — (a) An optical image of DM 5 from the 
Digitized Sky Survey (DSS). The data are taken 
from the SkyView web page. A white triangle 
denotes the position of the brightest far-infrared 
source (DR15 FIR 1). White circles represent 
the positions of catalogued infrared dark clouds, 
G79. 27+0.38 and G79. 34+0.33 . (6) A map of 
the CI 3 Pi- 3 Po emission integrated over the veloc- 
ity range Vlsr = —7.5 to +7.5 km s" 1 , smoothed 
with the FWHM=120" Gaussian weighting func- 
tion. Contours are set at intervals of 0.5 K 
km s _1 . (c) A map of CO J=3-2 emission in- 
tegrated over the velocity range Vlsr = —7.5 to 
+7.5 kms-\ smoothed with the FWHM=120" 
Gaussian weighting function. Contours are set at 
intervals of 1.5 K km s _1 . 

Fig. 2.— Observed spectra of CI 3 Pi- 3 P , CO 
J=3-2, 13 CO J=l-0, C 18 J=l-0 lines toward 
the CI peak (a), the centers of G79. 37+0.33 (6), 
G79. 27+0.38 (c), and the CO J=3-2 peak (d). 

Fig. 3. — (a) A map of the 13 CO J=l-0 emission 
integrated over the velocity range Vlsr — —7.5 
to +7.5 km s _1 , smoothed with the FWHM=45" 
Gaussian weighting function. Contours are set at 
intervals of 2 K km s _1 . White contours begin 
at 26 K km s" 1 . (b) A map of the C 18 J=l-0 
emission. Velocity range and smoothing width are 
the same as panel a) . Contours are set at intervals 
of 0.8 K km s _1 . White contours begin at 7.2 K 
km s" 1 . (c) Map of the 13 CO J=l-0 integrated 
intensity overlaid on the DSS image. Contours 
are set at intervals of 4 K km s _1 . (d) Map of the 
C 18 J=l-0 integrated intensity overlaid on the 
DSS image. Contours are set at intervals of 1.6 K 
km s . 

Fig. 4.— (a) Map of the CI 3 Pi- 3 P integrated in- 
tensity (contour) overlaid on the map of the 13 CO 
J=l-0 integrated intensity (gray scale), (b) Map 
of the [CI] integrated intensity (contour) overlaid 
on the map of the C 18 J=l-0 integrated inten- 
sity (gray scale). 

Fig. 5. — (a) Plot of the [CI] integrated intensity 
vs. the 13 CO J— 1-0 integrated intensity, (b) Plot 
of the [CI] integrated intensity vs. the C 18 J=l 
integrated intensity, (c) Plot of the C 18 J=l-0 
integrated intensity vs. the 13 CO J=l-0 J=l-0 
integrated intensity, (d) Plot of the [CI] integrated 



intensity vs. the CO J=3-2 integrated intensity. 

Fig. 6. — Comparison of the [CI] (filled cir- 
cles), CO J=3-2 (filled triangles), and C 18 
J=l (open triangles) integrated intensities. 
The cross section passes DR15A (Collcy 1980), 
(ai950,<W = (20 h 30 m 40.22 s ,+40°05'47.3"), 
and the [CI] peak position, (01950, <5i95o) = 
(20 h 30 m 36.0 s , +42°12'37.3"). 
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Table 1 
Telescope parameters 



Line 


Transition 


Frcq. [GHz] 


Telescope 


HPBW 


TVs [K] 


Vfss (^moon) 


Vmb 


Spectra 


CI 


3 Pi- 3 Po 


492.160651 


Mt. Fuji 


2'. 2 


900 


0.72 


0.45t 


539 


CO 


J=3-2 


345.795989 


Mt. Fuji 


3'.1 


500 


0.75 


0.67 1, 


324 


13 co 


J=l-0 


110.201353 


NRO 45 m 


17" 


400 


0.58 


0.48 


1088 


c 18 o 


J=l-0 


109.782182 


NRO 45 m 


17" 


400 


0.58 


0.48 


624 



tThc effective main-beam efficiencies including the effect of optical depth and the average temperature 
difference between the chopper inside the radome and the outside atmosphere. 



Table 2 
Column Densities 





T cx [K] 


N(CO) [cm" 2 ] 


N(C°) [cm- 2 ] 




iV(CI)/JV(CO) 


CI peak a 


18.0 


2.0xl0 18 


2.8xl0 17 


21-30 


0.15 


G79.37+0.33 b 


22.8 


8.4xl0 18 


1.9xl0 17 


87-130 


0.022 


G79.27+0.38 c 


18.1 


4.1 xlO 18 


7.2 xlO 16 


44-64 


0.017 


CO 3-2 peak d 


37.8 


3.3xl0 18 


9.3xl0 16 


51 


0.028 



a (ai950,<5i95o) = (20 h 30 m 36.0 s ,+40°12'37.3"). 
b (ai95oA95o) = (20 h 30 m 33.1 s ,+40°09'20.3"). 
c (ai950,<5i95o) = (20 h 30 m 07.9 s ,+40°08'12.3"). 
d (ai95o,<5i95o) = (20 h 30 m 43.9 s ,+40°06'37.3"). 

^Estimated from the CO column density using empirical relations, N(C 18 0) = 2.5 x 
10 14 (A V - 1.5) (Cernicharo & Guelin 1987) and iV(H) = 1.6 x 10 21 A V . 
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